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ABSTRACT 

Archival Hubble Space Telescope WFPC2 images were used to search for optical emission associated with the 
radio jets of a number of luminous quasars. From this search, we report new optical hotspot detections in the well- 
known blazar 3C 454.3 and the lobe-dominated quasars 3C 275.1 and 3C 336. We also find optical emission in the 
vicinity of the hotspot in 3C 208, but believe this is a chance alignment. Optical emission from the arcsecond-scale 
jet in 3C 454.3 is also detected. Multi-frequency archival radio data from the VLA and MERLIN are analyzed, 
and the synchrotron spectra of these high-power hotspots are presented. We estimate that their break frequencies 
are in the range of 10***- 10'^ Hz, with large uncertainties due to the wide gap in frequency coverage between the 
radio and optical bands. We also calculate their equipartition magnetic fields, and find that the anti-correlation 
between break frequency and magnetic field found by Brunetti et al. for lower power hotspots extends to these 
high power hotspots. This supports their model of hot-spots based on shock acceleration and synchrotron losses. 

Subject headings: Galaxies: active — galaxies: jets — quasars: general — quasars: individual (3C 208, 

3C 275.1, 3C 336, 3C 454.3) — radio continuum: galaxies 


1. INTRODUCTION 

The production of very high energy particles via diffusive 
shock acceleration has applications to many astrophysical prob¬ 
lems (e.g., Blandford and Eichler 1987). For instance, it is be¬ 
lieved to be responsible for the production of TeV energy par¬ 
ticles observed in supernova remnants, and in the very inner 
regions of active galactic nuclei (Catanese & Weekes 1999). 

Another important site of high energy particle acceleration is 
in the terminal hotspots of powerful extended extragalactic ra¬ 
dio sources. Their synchrotron spectra are observed from radio 
to optical frequencies (e.g., Meisenheimer et al. 1989), with 
the most energetic radiation requiring electrons with Lorentz 
factors, 7 ^ 10^, which must be accelerated in-situ. The re¬ 
cent rejuvenation of interest in the physics of hotspots is due 
mainly to the launch of the Chandra X-ray Observatory. Chan¬ 
dra’s sub-arcsecond resolution imaging capability and highly 
sensitive detectors allow for a unique probe of their physical 
conditions - e.g., magnetic field strength, when the X-ray emis¬ 
sion can be attributed to synchrotron self-Compton radiation, 
and energetics (e.g., Hardcastle et al. 2002). 

Prior to the launch of Chandra, most studies of hotspots con¬ 
centrated on high-resolution radio to optical mapping of their 
spectra, which can in principle, constrain their key physical pa¬ 
rameters - synchrotron break frequency, magnetic field, elec¬ 
tron mean free path. However, the higher frequency observa¬ 
tions were limited to ground-based optical and infrared instru¬ 
ments, and met with mixed success (e.g., Meisenheimer et al. 
1989, 1997). 

Recently, Brunetti et al. (2003) used the VET and success¬ 
fully detected optical counterparts to 6 out of 10 radio hotspot 
targets (one “non-detected” hotspot was limited mainly by con¬ 
fusion with a field source). Their targets were lower radio 
power hotspots than the more powerful ones which were pre¬ 
viously searched for optical emission (e.g., Meisenheimer et 
al. 1989, 1997). The lower power hotspots suffer less losses, 
so are sites of more efficient particle acceleration to the high 


energies ( 7 ~ 10 ^) which are required to produce optical emis¬ 
sion (Prieto et al. 2002). Conversely, the higher radio power 
hotspots, i.e. with larger magnetic fields, suffer greater (syn¬ 
chrotron and inverse Compton) losses, so show cutoffs in their 
spectra at lower frequencies - see the comparison between the 
spectra of the archetype high and low power hotspots in Cygnus 
A and 3C 445, respectively, in Figure 2 of Brunetti et al. (2003). 

In this scenario, there is a strong magnetic field dependence 
of the break frequency in their synchrotron spectra ( 17 , (xB~^) 
and observations of the lower power hotspots support this pic¬ 
ture (Brunetti et al. 2003). It is important to test this pro¬ 
posal, by extending these results to higher power hotspots (with 
greater B-fields). However, because their spectra are expected 
to cut off at lower energies, the higher power hotspots are faint 
in the optical, and can only be detected in deep Hubble Space 
Telescope (HST) images. 

The HST has demonstrated an unique capability of being 
able to detect faint optical counterparts to radio jet/hotspot fea¬ 
tures. This has proved to be important for the interpretation 
of the bright X-ray jets detected by Chandra (Tavecchio et al. 
2000; Celotti, Ghisellini, & Chiaberge 2001; Sambruna et al. 
2004). We were motivated by these work and after a success¬ 
ful detection of a bright optical knot in the well-studied quasar 
3C 279 from archival HST data (Cheung 2002; also motivated 
by the Chandra results on this object; Marshall et al. 2003), 
we initiated a similar search for extended optical features in 
other radio sources. Our search yielded optical detections of 
hotspots in three double-lobed 3CR quasars (3C 263, 3C 275.1, 
and 3C 336), the hotspot and jet knots in the blazar 3C 454.3, 
one tentative hotspot detection in 3C 208, and confirmation of 
the previously detected optical jet in PKS 1229-021 (Le Brun 
et al. 1997). The optical hotspot detection in 3C 263 has since 
been independently published by Hardcastle et al. (2002) who 
analyzed it in detail along with the Chandra detection. Neither 
3C 263 or PKS 1229-021 will be discussed further. In this pa¬ 
per, we focus our discussion on the results on the other hotspots 
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in light of the work of Brunetti et al. (2003). Initial results on 
3C 454.3 and PKS 1229-021 were presented in Cheung et al. 
(2003). 

In order to facilitate a comparison of our new optical detec¬ 
tions with those of Brunetti et al. (2003), we assume Hq = 50 
km s“' Mpc“^ and qo = 0.5 throughout this paper. 

2. SEARCH STRATEGY, DATA AND ANALYSIS 

The on-line Multimission Archive at the Space Telescope 
(MAST^) Science Institute was utilized to identify candidates 
for this study. MAST allowed us to identify archival HST im¬ 
ages of quasars by their radio flux by referencing the catalog 
of Veron-Cetty & Veron (1996). For this study, we chose an 
integrated 5 GHz flux greater than 0.5 Jy in order to ensure 
that the source contained radio features bright enough for con¬ 
vincing optical detections (assuming a typical range of radio- 
to-optical spectral indices of 0.8-1, where oc e.g., Sam- 
bruna et al. 2004). Our query in Summer 2002 resulted in about 
300 sources that had data obtained with the various imaging 
CCDs aboard HST. We focused on about 1/3 of these sources 
which had deep WFPC2 exposures (on both the high resolution 
planetary camera (PC) chip and any of the three lower resolu¬ 
tion wide field (WF) chips) using one or more orbits of total 
exposure time. It was in this type of exposure that the opti¬ 
cal jet in PKS 0637-752 was detected (Chartas et al. 2000; 
Schwartz et al. 2000). Lastly, we cross-referenced these candi¬ 
dates with known arcsecond-scale radio jets, aided by the cata¬ 
log of Liu & Zhang (2002) which led to about a dozen sources 
searched. This simple procedure was adopted in the interest of 
efficiently finding as many previously unidentified optical jets 
and hotspots as possible, with little intention of studying any 
complete sample. 

2.1. Archival Hubble Space Telescope WFPC2 Data 

Final combined and cosmic-ray removed HST WFPC2 im¬ 
ages were available from the WFPC2 Associations website (see 
footnote 3). Where there was no WFPC2 Associations image 
available, we downloaded the original exposures from the HST 
archive and combined them with the CRREJ task in IRAF^. The 
final data used for this study are listed in Table 1. 

Our identifications were made by overlaying high quality ra¬ 
dio maps (see below) in digital form with the optical data. HST 
has the well-known advantage over ground based telescopes of 
being very sensitive to faint point sources, so this procedure 
was essential to distinguish between field sources and radio jet 
related optical features. We also processed the optical images 
with an unsharped masking filter where smoothed versions of 
the original images were subtracted from the originals. This ad¬ 
ditional processing helped us to identify the nearest lying fea¬ 
tures to the bright optical nuclei (within ^1") in previous work 
with similar data (Cheung 2002; Sambruna et al. 2004), and led 
here to the successful detection of optical emission in the jet of 
3C 454.3 (§3.2). 

Because of this heightened sensitivity to point sources, many 
of the imaging fields were littered with sources unassociated 
with the radio source which could be mistakened for optical 
counterparts to jet knots. Therefore, in order to get a qualitative 
sense of the probability of a false detection, we made two tests. 
First, we repositioned the radio core onto several of the other 


prominent optical sources (field stars, galaxies) in the images 
and did not find any objects aligning with peaks in the radio 
jets within 1-2 WFPC2 pixels (i.e., to better than ^0.1" to 0.2"; 
each of our reported detections satisfy this). Second, we simply 
rotated the radio image in a random fashion with respect to the 
optical frame and again, did not find any convincing alignments 
with any knots in the radio jet. We therefore believe that the de¬ 
tections claimed here are real and are not chance alignments of 
field sources. 

Counts were measured from the optical hotspots, which were 
well-isolated from any other bright field sources, using circular 
apertures with the QPHOT task in IRAF. This gives essentially 
the counts in an infinite radius aperture by estimating the back¬ 
ground signal in circular annuli around the source. Counts in 
the features close to the bright nucleus in 3C 454.3 were mea¬ 
sured with square apertures. We estimated the contamination 
due to scattered light from the nuclei using identical apertures 
at other nearby positions on the CCD at equal distances from 
the core. These measurements suffered more contamination 
than in the hotspots so uncertainties were estimated accord¬ 
ingly by examining the fluctuations in the background. Count 
rates were converted to flux densities utilizing the inverse sen¬ 
sitivity measurements contained in the PHOTFLAM keyword 
in the image headers at the frequency given by the pivot wave¬ 
length indicated in the PHOTPLAM keyword. The WFPC2 fil¬ 
ters (Holtzman et al. 1995) closely approximate the standard 
optical bands: V (F606W), R (F622W, F675W, F702W), and I 
(F814W). The resultant measurements are reported in Tables 2, 
and 3, and extinction corrections were estimated using the near¬ 
est standard band values reported in the NED database derived 
from Schlegel, Einkbeiner, & Davis (1998). 

2.2. Archival Radio Data 

The quasars in which we successfully identified optical emis¬ 
sion in the jets and/or hotspots are well-known, and have been 
extensively imaged at centimeter wavelengths, so this analy¬ 
sis did not warrant new observations. High quality radio ob¬ 
servations were obtained from the archival databases of the 
NRAO Very Large Array (Thompson et al. 1980), and ad¬ 
ditional MERLIN observations for the three lobe-dominated 
quasars. The MERLIN data were dual-frequency synthesis L- 
band observations which utilized a seven-element array of Tab- 
ley, MK2, Darnhall, Wardle, Knockin, Cambridge, and Def- 
ford. 

The list of archival radio data collected are summarized in 
Table 1 with some effort made to find the original references. 
We invariably selected 5 GHz A-configuration data (0.4" to 
0.5" resolution) for each source. The radio observations are all 
to our knowledge, unpublished. In one instance (3C 454.3), we 
used 8.4 and 22 GHz data from one of our own recent observ¬ 
ing programs when the source was used as a phase calibrator for 
two targets. The VLA data were calibrated using standard pro¬ 
cedures in the NRAO AIPS reduction package (Greisen 1988; 
Bridle & Greisen 1994) and exported to the Caltech DIEMAP 
program (Shepherd, Pearson, & Taylor 1994) for imaging and 
deconvolution. The initial calibration of the MERLIN data 
were performed at Jodrell Bank by the archivists and then self- 
calibrated and imaged at Brandeis. The 1.4 and 1.7 GHz MER¬ 
LIN datasets were imaged and analyzed separately and the final 
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Summary of Multi-Telescope Archival Data 


Source 

(1) 

Instrument 

(2) 

Date 

(3) 

Frequency 

(4) 

Exp Time 
(5) 

Program 

(6) 

Observer or Reference 

(7) 

3C 208 

MERLIN 

14 Aug 1993 

1.413+1.658 

14 hr 

- 

- 


VLA A-configuration 

13 Dec 1992 

8.44 

3500 

AL280 

S. Garrington 


HST WFPC2 WF3/F675W 

23 Mar 1997 

4.463x10^ 

9400 

6491 

A. Stockton 

3C 275.1 

MERLIN 

07 Aug 1993 

1.413+1.658 

12 hr 

_ 

_ 


VLA A-configuration 

06 Mar 1982 

4.89 

2340 

WARD 

J. F. C. Wardle 


VLA A-configuration 

13 Dec 1992 

8.44 

3490 

AL280 

S. Garrington 


ffST WEPC2 PC1/E675W 

25 Jul 1995 

4.463 xlO^ 

1800 

5978 

S. Rawlings 

3C336 

MERLIN 

13 Aug 1993 

1.413+1.658 

15 hr 

_ 

_ 


VLA A-configuration 

06 Mar 1982 

4.89 

2310 

WARD 

J. F. C. Wardle 


VLA A-configuration 

13 Dec 1992 

8.44 

4440 

AL280 

S. Garrington 


HSTWFPC2 WF2/F814W 

20 Jul 2000 

3.755 xlO^ 

6600 

8220 

M. Franx 


HST WFPC2 WF2/F702W 

12-14 Sep 1994 

4.334 xlO^ 

24000 

5304 

Steidel et al. (1997) 


HST WFPC2 WF3/F622W 

20 Dec 1994 

4.846x10^ 

7200 

5401 

Ridgway & Stockton (1997) 


HST WFPC2 WF2/F606W 

20 Jul 2000 

4.993 xlO^ 

6600 

8220 

M. Franx 

3C 454.3 

VLA A-configuration 

31 Jan 1985 

1.50 

10860 

AC 120 

R. Perley 


VLA A-configuration 

31 Jan 1985 

4.86 

10490 

AC 120 

R. Perley 


VLA B-configuration 

06 Jun 2002 

8.46 

860 

AC641B 

C. C. Cheung 


VLA BnA-configuration 

28 May 2002 

22.46 

3850 

AC641A 

C. C. Cheung 


HST WFPC2 WF2/F702W 

06 Dec 1997 

4.333 xlO^ 

2100 

6619 

Chen et al. (2001) 


Table 1 

Column (1) specifies the common name of the sources. 

Column (2) indicates the telescope and instrument/configuration used. The HST WFPC2 CCD chip and filters used are indicated. 

Column (3) is the UT date of the archival observation. 

Column (4) is the observing frequency in GHz. 

Column (5) is the time on source in seconds unless otherwise indicated. The integration time for the radio data correspond to the unedited data. 

Column (6) are the VLA or HST program codes/numbers. 

Column (7) is reference for observation or observer name if we could not locate the reference. 
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Lobe-dominated 3CR Quasar Hotspot Multi-wavelength Flux Densities 


Frequency 

Flux 

Flux 

Flux 

(GHz) 

3C 208-West 

3C 275.1-North 

3C 336-South 


1.53 

0.108±0.016 

0.355±0.053 

0.097±0.015 

4.89 

0.0277±0.003‘' 

0.178±0.027 

0.045±0.007 

8.44 

0.021 ±0.003 

0.106±0.016 

0.029±0.004 

3.755 xlO^ 



0.18(0.2 l)±0.04xl0-<’ 

4.334x10^ 

4.463x10^ 

0.21*(0.23)±0.04xl0-® 

0.30(0.32)±0.06xl0-‘’ 

0.15(0.17)±0.03xl0“<’ 

4.846x10^ 



0.11(0.13)±0.03xl0-<’ 

4.993x10^ 



0.14(0.17)±0.05xl0-<’ 


^radio 

1.02±0.14* 

0.69±0.12 

0.74±0.12 

Otro 

1.03 (1.02)* 

1.16(1.16) 

1.11 (1.10) 


Table 2 

Fluxes are in units oe Janskys with extinction corrected values indicated in parentheses. 

Radio spectral indices, (F oc are best-eit linear fits to the radio data. The radio-to-optical spectra, a„ are 
CALCULATED USING THE 5 GHZ AND 4.463 X 10^ GHZ MEASUREMENTS FOR 3C 208 AND 3C 275.1, AND 4.334 X 10^ GHZ FOR 3C 336. 
A SINGLE POWER LAW FIT THROUGH ALL OF THE DATA POINTS YIELD SPECTRAL INDICES OF 1.04, 1.14, AND 1.08, FOR THE THREE 

SOURCES, RESPECTIVELY. 

“ Taken from Bridle et al. (1994). 

The value indicated with an asterisk (*) is an optical candidate. 


3C 454.3 Multi-wavelength Flux Densities 

Frequency 

Flux 

Flux 

Flux 

Flux 

(GHz) 

Knot A 

Knot B 

Knot C 

Hotspot 

1.50 




0.541±0.054 

4.86 




0.222±0.022 

8.46 

0.013 ±0.004 

0.008±0.002 

0.010±0.003 

0.132±0.013 

22.46 




0.065±0.010 

4.333 xlO^ 

0.25(0.32)±0.25xl0 

0.25(0.32)±0.13xl0-® 

0.12(0.15)±0.06xl0-^ 

0.50(0.63)±0.10xl0-‘’ 

^radio 




0.80±0.06 

CXfo 

1.00 (0.98) 

0.96 (0.93) 

1.04(1.02) 

1.15(1.13) 


Table 3 

Fluxes are in units of Janskys with extinction corrected values indicated in parentheses. 

The radio spectral index for the hotspot, (F « is best-fit linear fits to the radio data. The radio-to-optical 
SPECTRA, Oro ARE CALCULATED USING THE 8.4 GFIZ AND 4.333 X 10^ GFIZ MEASUREMENTS. 
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images averaged. 

We used DIFMAP’s MODELFIT program to model the 
brightness distribution in the hotspot regions in the radio 
sources with elliptical gaussian components. The inner jet in 
3C 454.3 was only detected and resolved adequately in the 8.4 
GHz data so radio fluxes were measured for the jet only at 
this frequency. The hotspot regions in the lobe-dominated 3CR 
quasars are complex and our estimate of the errors in the ex¬ 
tracted fluxes are about 15%. The 5 GHz flux for the 3C 208 
hotspot was taken from Bridle et al. (1994). 

The resultant fluxes extracted in this manner are presented 
in Tables 2 and 3. The radio spectral indices (ar) are linear 
best fits to the data and are presented along with the radio-to- 
optical spectra indices (aro)- In 3C 208, the radio spectrum of 
the hotspot is notably steeper than the other three cases where 
we claim bona-fide optical detections. The optical flux of the 
optical source near the radio hotspot lies on a power-law ex¬ 
trapolation of the radio spectrum (i.e., ~ 0 ;^) so it may still 

be due to synchrotron radiation if the spectrum breaks at even 
higher energies. In the other three cases, the radio-to-optical 
spectra are greater than the radio only spectra, indicating steep¬ 
ening of the hotspot spectra in the observed frequency range. 

3. RESULTS AND DESCRIPTIONS OE INDIVIDUAL SOURCES 

We successfully identified previously unknown optical coun¬ 
terparts to the hotspots of the lobe-dominated quasars 3C 275.1 
(§3.1.2; Figure 2) and 3C 336 (§3.1.3; Figure 3), and in the 
blazar 3C 454.3 (§3.2; Figure 4). We also identified a possible 
optical counterpart to the hotspot in the lobe-dominated quasar 
3C 208 (§3.1.1; Figure 1). As mentioned in the introduction, 
the optical hotspot detection in 3C 263 (z=0.652) was reported 
since by Hardcastle et al. (2002) along with the new Chandra 
X-ray detection so is not discussed here. The lobe-dominated 
quasars are best known to many as being in the Hough & Read- 
head (1989) sample of double lobed 3CR quasars. A recent 
summary of their source properties and VLBI observations is 
in Hough et al. (2002). 

3.1. The Double-lobed SCR Quasars 

3.1.1. Optical Emission Nearby the Radio Hotspot in SC 208 

The optical identification of the western radio hotspot in 
3C 208 (z=l.ll) is only tentative since the optical feature is 
not as precisely-aligned with any peak in the radio structure 
(Figure 1) as in the other cases presented in this paper. The 
closest radio feature is knot “B” in the nomenclature of Bri¬ 
dle et al. (1994) which we identify as the hotspot. The optical 
feature appears resolved; although the radio and optical peaks 
are ~0.3±0.1" apart (1.6 kpc projected), some of the optical 
light may still be associated with the radio hotspot. We could 
not confirm the optical detection in shallower WFPC2 F555W 
and F702W images (snapshots) taken on the high resolution 
PC chip (Lehnert et al. 1999). There is no apparent optical 
emission from the radio jet in any of the images. If all of the 
optical emission near the 3C 208 radio hotspot is indeed asso¬ 
ciated with the radio source, its overall spectrum is consistent 
with a single power-law with little or no steepening from the 
observed radio spectral index of ~1. Further observations with 
multiple optical filters to measure the spectrum of this feature 
are encouraged. 

3.1.2. Hotspot Detection in SC 275.1 


3C 275.1 (z=0.557) is a radio source with a “dog-leg” ap¬ 
pearance (Stocke, Burns, & Christiansen 1985). A MERLIN 
1.4 GHz image appears in Akujor et al. (1994) and a deep 
VLA 5 GHz image in Stocke et al. (1985). The optical field 
in 3C 275.1 has been well studied (e.g., Hintzen & Romanishin 
1986, and references therein). 

Optical emission from the northern hotspot is prominent in 
the unprocessed HST image (Eigure 2) and is clearly coincident 
with the peak in the high resolution MERLIN map. This detec¬ 
tion was independently confirmed in a recent study by Hardcas¬ 
tle et al. (2004). It was never noted in previous ground based 
optical imaging studies presumably due to its faintness. There 
is some extended optical emission around the hotspot, some 
of which is aligned along the complex radio jet/hotspot struc¬ 
ture. The hotspot has also been detected in the X-rays with 
Chandra at a flux level above the extrapolation of the radio-to- 
optical spectrum signaling an additional spectral component, 
presumably inverse Compton emission (Crawford & Eabian 
2003; Hardcastle et al. 2004). 

3.1.3. Hotspot Detection in SC SS6 

3C 336 (z=0.927) is a large angular size source, stretching 
across almost 1/2 arcmin in the sky (~250 kpc projected). It 
lies in an extremely dense optical field (Steidel et al. 1997). In 
the nomenclature of Bridle et al. (1994), the optically detected 
feature (Eigure 3) is knot “C” in the jet to the south. Although 
it marks the termination point of the jet, the feature does not 
conform to their strict classification as a hotspot. Rather, they 
classify a radio feature further south as the hotspot although C 
appears more compact - we will proceed with our discussion 
assuming that C is indeed a hotspot. 

Bridle et al. (1994) found the 5 GHz flux to be 38.3 mJy 
compared to our observed 45 mJy. Our fluxes were not cor¬ 
rected for contribution of the surrounding lobe emission so may 
be systematically higher than the fluxes derived by Bridle et al. 
(1994). Region C is actually resolved into a peak where the 
optical feature is clearly coincident, and a weaker extension to 
the west. Our flux density measurement is from the dominant 
feature. Eainter extended optical emission may be associated 
with the radio tail but it is at the same level as surrounding low 
level emission which is clearly not associated with the radio 
structure. 

We found the optical hotspot in four separate WEPC2 im¬ 
ages obtained in different programs (Table 1). The published 
figure of the deep 24,000 second WEPC2 P702W exposure in 
Steidel et al. (1997) actually shows the feature (at the bottom 
edge of their figure 2 d), but was not identified by those authors 
as associated with the radio jet. This is the only hotspot where 
we could measure an optical spectrum, aopucai ~ 1.5 ±1.0, but 
the spectral coverage in the optical is limited making this value 
highly uncertain. The difference between the radio and optical 
indices is formally, Aa ~0.8. 

3.2. The Optical Hotspot and Jet in SC 454.S 

3C 454.3 (z=0.859) is a superluminal 7 -ray blazar with a 
core-jet VLBI morphology (Pauliny-Toth et al. 1987; Jorstad 
et al. 2001). The parsec-scale jet stretches westward from the 
nucleus, continuing with more extended emission to the north¬ 
west. This extended emission traces a path toward the kilo- 
parsec scale structure seen in VLA images, where the opti¬ 
cal emission was detected in the HST data (Eigure 4). A re¬ 
cent Chandra image detected X-ray emission from the jet and 
hotspot (Marshall et al. 2005). 



6 


Cheung et al. 



dUIn FLA. fn 


aCSDS Hotpot 

H 



H -n 

RA [«t44] 


Fig. 1.— Radio (contour) and optical F675W (colorscale) overlay of 3C 208 with large field of view (top) and zoomed in view centered near the western hotspot 
region (bottom panel). The top image shows the VLA 8.4 GHz image convolved to 0.25^^ beam. The base level is 35/iJy/beam and spaced by factors of 2. The 
nucleus is marked with a white cross. The MERLIN 1.4 GHz image on the bottom panel has a 0.231^' x0.16F^ beam at PA=22° (uniform weighting). The minimum 
contour is 1.5 mJy/beam spaced by factors of 2. Note that the radio/optical features in the hotspot region are not quite aligned making this only a tentative optical 
hotspot identification. North is up and east is left on the figures. 
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Fig. 2.— Radio (contour) and optical F675W (colorscale) overlay of 3C 275.1 with large field of view (left) and zoomed in view centered near the northern 
hotspot region (right panel). The left image shows the VLA 8.4 GFIz image convolved to 0.25” beam. The base level is 0.125 mjy/beam and spaced by factors of 2 
up to seven contours to show the optical hotspot emission more clearly. The nucleus is marked with a white cross. The MERLIN 1.4 GHz image on the right panel 
has a 0.2''x0.15” beam at PA=34° (uniform weighting). The minimum contour is 0.75 mJy/beam spaced by factors of 2. North is up and east is left on the figures. 


A deep WFPC2 image of this quasar was obtained as part 
of a program to search for Ly-a absorption systems in its field 
(Chen et al. 2001), so the exposure was centered on one of 
the low resolution wide-held (WF) CCD chips (about twice the 
pixel size of the higher resolution PC chip). An unresolved 
(FWHM~0.15-0.2") optical source coincident with the posi¬ 
tion of the hotspot in the radio jet is apparent in the HST image, 
very near one of the diffraction spikes from a nearby bright 
point source (Figure 4; see also the unsharp-mask processed 
image in Figure 1 of Cheung et al. 2003). 

Upon closer inspection, a series of three bright optical knots 
are present which match peaks in the jet as seen in the full res¬ 
olution 8.4 GHz radio map (a hybrid BnA conhguration dataset 
with extra north-south resolution). These radio features are 
conhrmed in a MERLIN 1.7 GHz map published in Browne 
et al. (1982). Our preliminary processing (Cheung et al. 2003) 
of the HST data used a gaussian smoothing with a very broad 
profile so that a large negative bowl suppressed emission from 
these inner knots. 

The optical feature coincident with radio knot B is seen in a 
280s WFPC2 PC snapshot image obtained as part of the Lehnert 
et al. (1999) survey of the 3CR sample. The feature is clearly 
separated from the nucleus and both its position and flux (in the 
same hlter as in the current deep WF image) almost match ex¬ 
actly. The knot is about a 3cr detection. The fainter outer knot 
(C) could not however be verihed in this image as expected 
from signal-to-noise limits, and knot A is simply too close to 
the bright nucleus to be distinguished easily from the scattered 

^ See: http://www.stsci.edu/instraments/wfpc2/ 


light. Although the integrated flux of the hotspot is greater than 
knot B, it was undetected in the higher resolution (PC chip) 
snapshot image. For the measured peak surface brightness from 
the WF image of about 6 /iJy/arcsec^, we hnd that the expected 
SNR is at best about unity for the PC chip exposure using the 
on-line WFPC2 exposure time calculator^. 

The jet in the 22 GHz data was very ill-defined so fluxes 
for the inner jet knots were not extracted from this dataset. 
However, the hotspot was prominent. Our total 22 GHz flux 
for 3C 454.3 was renormalized to contemporaneous single dish 
measurements from the Metsahovi quasar monitoring program 
(H. Terasranta, 2003, private communication; Terasranta et al. 
2004). The radio source varied by at most 3% over the 1 month 
period, so the average of five measurements obtained by the 
Metsahovi group was used. 

4. DISCUSSION 

Brunetti et al. (2003) recently found that lower power ra¬ 
dio hotspots are brighter at optical wavelengths than these more 
powerful radio hotspots which are customarily searched for op¬ 
tical emission (Meisenheimer et al. 1989). They propose a 
simple in-situ shock acceleration scheme whereby high radio 
power hotspots with larger B-flelds suffer greater losses, and 
are inefficient at accelerating electrons up to 7 ~ 10 ^, required 
to produce optical emission. In this scenario, the position of 
the “break” frequency in the observed spectrum will depend on 
the balance between the rate of electron acceleration in a shock 
region and combined synchrotron and inverse Compton losses. 
There is a strong magnetic field dependence of the break fre- 
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Fig. 3.— Radio (contour) and optical F702W (colorscale) overlay of 3C 336 with large field of view (left) and zoomed in view centered nearby the southern 
hotspot region (right panel). The left image shows the VLA 8.4 GHz image convolved with a 0.25'' beam. The base level is 0.15 mJy/beam and spaced by factors 
of 2 up to five contours to show the optical hotspot emission more clearly. The nucleus is marked with a white cross. The MERLIN 1.6 GHz image on the right 
panel has a 0.221" x 0.152" beam at PA=16.8° (uniform weighting). The minimum contour is 1 mJy/beam spaced by factors of 2. North is up and east is left on the 
figures. 


Properties of the Three 3CR Quasar Hotspots 


Source 

(1) 

R (asec, cm) 

(2) 

log(i4,) 

(3) 

Beq 

(4) 

3C 275.1-N 

0.22, 5.8 

10 .0±1.2 

2.6 

3C 336-S 

0.25, 5.9 

11 .1±1.2 

2.3 

3C 454.3-NW 

0.25, 6.4 

11.4±0.6 

3.8 


Table 4 

(2) The assumed radius of a spheroid in arcseonds ("), and in units oe 10^' cm. 

(4) The logarithm oe the derived break erequency, Ub, in Hz (see §4). 

(3) Equipartition magnetic field, Beq in units of 10^ G. These assume a uniformly filled region, negligible beaming, 

AND 'Ymin— 10. 
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Fig. 4.— Radio (contour) and optical HST F702W (colorscale) images of 3C 454.3 overlayed, and a zoomed in view of the one-sided jet (right). The nucleus is 
marked with a white cross. The radio image is from a VLA 8.4 GHz dataset convolved with a beamsize of 0.7" x 0.3'' at a PA=67° (uniform weighting). The base 
level is 1 mJy/beam and spaced by factors of 2; only six contours are plotted on the left panel to show the optical hotspot emission more clearly. The HST image 
is shown before (left) and after (right) the subtraction of the elliptical isophote model of the point source. The four diffraction spikes around the optical nucleus are 
clearly present, as is a portion of a spike in the upper right of the displayed field of view from a field star (left panel). North is up and east is left on the figures. 


quency in their synchrotron spectra, i^h (Figure 6), and 

observations of the lower power hotspots support this predic¬ 
tion (Brunetti et al. 2003; see also earlier discussions of this 
scenario by Meisenheimer et al. 1989 and Prieto et al. 2002). 

It is important to extend this work to higher power hotspots 
(with greater B) which are expected to be faint in the near-IR 
and optical, and can only be detectable in deep HST expo¬ 
sures. Our detections of the faint optical counterparts in the 
three sources help us to constrain in the spectra of these high 
power hotspots. 

Figure 5 shows the radio-to-optical spectral energy distribu¬ 
tions (SED) of the three hotspots with bona-fide optical de¬ 
tections from our analysis of the archival HST and radio data. 
The optical fluxes clearly lie below a simple extrapolation of a 
power law from the observed radio data, so are consistent with 
steepening in the high energy spectrum. 

The data are admittedly sparse and the overall shape of the 
SED can only be roughly constrained. We assume that the 
SED can be approximated by a simple smoothly varying double 
power law function so that we can parameterize it with a break 
frequency (r^) at which the power law slope changes. The low 
energy slopes are set by the observed radio spectra (ar). In the 
absence of precise information on the high energy slopes, we 
assume that the optical spectral index, ao=ar-0.5, which is ex¬ 
pected from first order Eermi shock acceleration models with 
synchrotron losses (e.g. Heavens & Meisenheimer 1987). This 
is indeed, very near what is observed in the case of the 3C 336 
hotspot, where we have a rough constraint on the optical spec¬ 
trum (§3.1.3). In this manner, is related to observed variables 
by: log (i/b) = log (t'o) -i- 2 (aro-ai) log (t'r/r'o), where the sub¬ 


scripts indicate the radio (5 GHz) and optical values. In prac¬ 
tice, the uncertainty in determining is dominated by a{ai), 
since (T(ar) ~0.1, compared to cr(Q:ro)~0.02-0.03. Therefore, 
(T(log (i%)) 2 log (h'r/i'o) cr(o:r) ~ 10 o-(ar). Our estimates of 

Vh ~ 10'°- 10" Hz for these hotspots are therefore good to ± 
1 dex (Table 4). 

Eollowing Brunetti et al. (2003), we used the observed ra¬ 
dio spectra, flux densities at 5 GHz, and constraints on the 
sizes from our modelflts, to calculate magnetic fields assuming 
equipartition (Table 4). These magnetic fields are in the range 
of a few hundred ^G, so are relatively large. As in Brunetti 
et al. (2003), we have assumed for simplicity that the hotspot 
emission is not relativistically beamed (but see Georganopoulos 
& Kazanas 2003). If the radiation is beamed, the equipartition 
magnetic fields calculated must be reduced by a factor of 6^^’’ 
to transform to the source rest-frame (Stawarz, Sikora, & Os- 
trowski 2003), where 6 is the Doppler factor. 

Recent work by Sambruna et al. (2004) and Hardcastle et 
al. (2004) provide additional optical hotspot source detections 
and we were able to estimate (i^, B^q) for 8 additional hotspots 
in 6 radio sources as described above (Table 5). Although in¬ 
formation for a large number of optical hotspots are provided 
by Hardcastle et al. (2004), only four of these hotspots (in three 
objects) were not previously considered in the study by Brunetti 
et al. (2003). Of these four hotspots, we could not confirm the 
case in 3C 280 from our independent analysis of the same HST 
data (see also the radio/optical image overlay in Eigure 3 of 
Ridgway et al. 2004 - the radio and optical peaks were not 
clearly coincident as in our case in 3C 208; § 3.1.1) - this case 
is omitted in our discussion. 
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Optical Hotspots from the Literature 


Source 

(1) 

z 

(2) 

Gr 

(3) 

Ref. 

(4) 

CKro 

(5) 

Ref. 

(6) 

log(l4,) 

(7) 

Beq 

(8) 

0405-I23-N 

0.574 

0.90±0.10 

C04; S04; T05 

1.10 

C04; S04; T05 

12.7±1.0 

2.9 

3C 228-S 

0.5524 

0.76* 

J95 

1.05 

H04 

11 .8±1.0 

1.7 

3C 245-W 

1.029 

0.87±0.08 

C04; S04; T05 

1.39 

C04; S04; T05 

9.5±0.8 

6.0 

1136-135-W 

0.554 

0.85±0.08 

C04; S04; T05 

1.22 

C04; S04; T05 

11 .0±0.8 

3.3 

1150-F497-H 

0.334 

0.72±0.09 

C04; S04; T05 

1.12 

C04; S04; T05 

10.7±0.9 

1.2 

I150-F497-I 

0.334 

0.80±0.10 

C04; S04; T05 

1.07 

C04; S04; T05 

12 .0±1.0 

1.4 

3C 403-PI 

0.059 

0 .8* 

D99 

0.87 

H04 

14.0±1.0 

1.1 

3C 403-P6 

0.059 

0 .8* 

D99 

0.85 

H04 

14.1±1.0 

1.4 


Table 5 

‘Errors of 0.1 assumed. 

(3) Radio spectral index from references in column (4). 

(5) Radio-to-optical spectral index from references in column (6). We used fluxes reported by the authors at 5 GHz 
IN THE RADIO (IN 3C 228 AND 3C 403, H04 CONVERTED FROM 8.4 GHZ USING ar=0.5), AND 3.45 X lO*'* (3C 228), 4.34x lo'"^ 

(3C 403), AND 5.13x10''* Hz (rest) for the optical. 

(7) The logarithmic of the derived break frequency, Ub, in Hz (see §4). 

(8) Equipartition magnetic field, Beq IN UNITS OF 10^ G. See Table 4 for assumptions used. 

References in columns (4) and (6): C04 = Cheung (2004), Cheung et al. in prep.; D99 = Dennett-Thorpeet al. (1999); 
H04 = Hardcastleet al. (2004); J95 = JOHNSON ET AL. (1995); L92 = LiU et al. (1992); S04 = Sambruna ET AL. (2004); T05 = 

Tavecchioet al. (2005) 


We find that the derived {v/,, B) of the high radio power 
hotspots (Tables 4 & 5) lie near the extrapolation of the Vh (xB~^ 
trend traced by Brunetti et al. (2003), which lends support 
to their simple picture (Figure 6). This is a remarkable find¬ 
ing considering the hotspots span about 3 orders of magnitude 
in magnetic field energy density, while considering both radio 
galaxies and quasars. In the scenario outlined by Brunetti et al. 
(2003), the proportionality is: i/b oc (m; where r 
is the dynamical age of the hotspot and the ratio of the energy 
densities of the lobe (m;) and the hotspot (mhs) accounts for adi¬ 
abatic losses. More precise mapping of the high-power hotspot 
SEDs are necessary to confirm that the -3 slope indeed extends 
over 6 orders of magnitude in frequency (e.g. that there is no 
change in the slope). If confirmed, these observations could im¬ 
ply that that there is a single dynamical age for radio hotspots 
in radio galaxies and quasars over a wide range in power. 

Spitzer Space Telescope observations, already approved for 
two of our new optical hotspot sources, will allow us to measure 
the high energy slope of the SED which will constrain i/h more 
robustly. While additional millimeter and sub-mm observations 
(in the future with ALMA) are required to map the shape of the 
high energy synchrotron spectrum, existing radio, optical, and 
new infrared data can constrain the break frequency adequately 
(to better than a decade) for the purpose of showing the extrap¬ 
olation of the Z 4 , - Beq sequence determined from the theory. 
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Eig. 5.— Spectral energy distributions of the high power hotspots detected in the two double-lobed quasars and the blazar 3C 454.3. Break frequencies ( 14 ,) are 
only roughly constrained with the limited observations by fitting a double power law spectrum to the data. 
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Fig. 6. — Break frequency (i^,) in the synchrotron spectra of the hotspots with optical detections vs. magnetic field calculated assuming equipartition (Beq); 
adapted from Brunetti et al. (2003). Our new optical hotspot detections in 3 powerful quasars (squares), and other detections from Hardcastle et al. (2004) and 
Sambruna et al. (2004) (triangles) follow approximately the i/b trend (solid line) expected from simple in-situ particle acceleration theory (Brunetti et al. 

2003). Uncertainties in L'b for the hotspots we analyzed (Tables 4 & ^ are indicated with error bars. For the rest, the break-frequencies are within a factor of 2 as 
reported by Brunetti et al. (2003), with the exception of the case of 3C 35IL where the line connecting two points indicates the uncertainty. 




